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1. Introduction 

The current stage of theoretical cosmology and, in particularly, the dark energy problem 
calls to consideration of new field theory models which maybe useful for the explanation 
of the universe acceleration. In the situation when the attempts to use the well- 
known classes of field theories as candidates for dark energy (dark matter) are not 
quite successful, such strategy seems to be justified. Moreover, one should bear in mind 
that number of problems of the early time cosmology (say, the fundamental cosmological 
constant problem) are still too far from the resolution. The (old or new) field theoretical 
model which should include also gravitation itself should solve these problems at once 
(at least, in the regime where one can still neglect the quantum gravity effects). 

In the recent paper pP the class of dark energy models where dark energy (scalar or 
spinor field theory, etc) couples with some power of the curvature has been proposed. It 
has been demonstrated on the example of scalar Lagrangian coupled with the power of 
curvature that such modification of gravity-matter coupling term may explain not only 
the current cosmic speed-up, but also the dark energy dominance assisted by gravity. On 
the same time, theories of similar nature have been recently proposed as good candidate 
for dynamical resolution of the cosmological constant problem E| (see also |3j). Of 
course, such models are not standard ones, in the sense that they are not multiplicatively 
renormalizable in curved spacetime Hence, they should be considered as kind of 
effective theories (without clear understanding of their origin and their relation with 
more fundamental string/M-theory). 

In the present paper we continue the study of the properties of such models on the 
example of scalar self-interacting theory non-linearly coupled with some power of the 
curvature. Admitting the discrete symmetry which prohibits odd powers of scalar in 
the potential, the effective action formulation (section 2) is developed in close analogy 
with multiplicatively renormalizable theories [5J. Using Riemann normal coordinate 
expansion, the modified Green function is evaluated at weak curvature limit (section 
3). (We work in linear curvature approximation what is expected to be appropriate for 
the study of quantum effects in the inflationary universe and in late-time, dark energy 
universe). Subsequently, this Green function is used for the calculation of the one-loop 
effective action in the same limit (section 4). After the one-loop renormalization of 
the effective action done in section 5, one arrives at finite expression which is used for 
the numerical study of the phase structure in FRW universe. The number of examples 
explicitly demonstrates the behaviour of the effective potential for such effective theory. 
In particularly, it is shown the possibility of curvature-induced phase transitions in the 
model under discussion (section 6). It is indicated that radiative corrections are not 
very relevant at weak curvature. 

Finally, in section 7 it is shown that scalar theory under consideration may serve 
also for dynamical resolution of cosmological constant problem. The explicit choice 
of parameters (which maybe considered as kind of fine-tuning) giving such dynamical 
resolution is presented. It is interesting that radiative corrections are not relative for 
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such a mechanism. Some discussion and outlook are given in the last section. 



2. 4 theory coupled non-linearly with curvature power 

We extend the 4 theory which is the simplest model where the spontaneous symmetry 
breaking takes place. One of the possible couplings of scalar Lagrangian with curvature 
(at some power) maybe included into the starting action: 



S = I d D x y/^g 



(i) 



2k 2 \M 2 

where M is an arbitrary mass scale, g is the determinant of the metric tensor and 
Ld is the ordinary Lagrangian density of the 4 theory, 



cm = - + f *g - ^5, (2) 

where 0o is a real scalar field. Our sign convention for the metric is (H • • •) and 

follows the notation of Birrell and Davies jH]- If the spacetime curvature is constant, 
the non-linear curvature coupling disappears by the transformation 



-a/2 



Ao-^VV (4) 



\ M2 , 

We are interested in a non-trivial effect of the non-linear curvature coupling. It will 
be found in a non-static and/or inhomogeneous spacetime. Note that the non-linear 
coupling of above sort has been proposed in Ref. [T] (see also related model in [2j) as 
the model of dark energy which naturally resolves the problem of dark energy dominance 
in the current universe. 

The Lagrangian is invariant under the discrete transformation, 

- -0- (5) 

This symmetry prevents the Lagrangian from having <ft 3 terms. A non-vanishing 
expectation value for the field breaks the discrete symmetry spontaneously. To 
evaluate the expectation value for we calculate the effective action. 

We start from the generating functional VK[J]. The generating functional of the 
theory is given by 

e mA= [v<f>oel s+ y dD *^*o J (6) 



In the presence of the source J the classical equation of motion becomes 
5S 
5(p (x 

In curved spacetime the functional derivative is defined by 



= ~J(x). (7) 

<t>0=4>b 



5[0 O + 50] - S[<f> ] = I d D x^^- ) 5Mx). (8) 
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We divide the field O m t° a classical background (fit, which satisfies the Eq.(0) and 
a quantum fluctuation 0, 

0o = 0b + fr 1/2 0. (9) 
In terms of <f> b and the action S[4>q] is rewritten as 

s[M = s[(f> b + h l ' 2 $\ 



S[(j) b ) + I d D x^~g 



55[0 C 



S(j) (x) 



h 1/2 <p(x) 



4>o=4>b 



+ 



h 



d D x J d D y^^^^)^x)tGf(x,y)4>(y) + 0(4> 3 ) 
d D x^=gJ(x)h l/2 4)(x) 



2 

+ ~ I d D x I d D y v / -g(x) v / -g(y)(p(x)iG d l (x,y)4)(y) 



+ O(^ 2 3 ), 
where G d 1 (x, y) is given by 

5 2 S[0 O ] 



(10) 



iG d x {x,y) = 



50 O (x)50o(y) 



R 



R 



M 2 



a-l 



R; 



M 2 



Therefore the generating functional (JHJ) is expanded to be 

X j X?0 e ll d D x J d D y^^g(x)^-g(y)4>(x)iG^ 1 (x,y)4>(y)+0(^^) ^) 

Performing the path-integral, 



we obtain the generating functional 



(13) 



ih 



W[J] = S[<f> b ] - J d D x^J{x)(t) b {x) - — In [Det iG d l ] + 0(fr 3/2 ). (14) 

The effective action r[0 c ] is given by the Legendre transform of VF[J], 

T[0 C ] = W[J] - J d D x^<f> c (x)J(x), (15) 

where C denotes the expectation value of O in the presence of the source J, 
, _ 5W[J] 



5J 



(16) 



Substituting Eq. (jT4*)l to Eq. (JT5|) we find that the effective action is expanded to be 



T[0 C ] = 5[0 C ] - y In [Det iG, 1 ] + 0(fr 3 / 2 ). 



(17) 



The one-loop effective action in 4> A theory coupled non-linearly with curvature 



5 



In the path integral formalism the expectation value of 0o hi the ground state 
satisfies 



SW[J] 

This equation is rewritten as 



o. 



0<Pc ^-,(00) 

where we have used the relation 

sn<t>c] _ j 



;is) 



(19) 



(20) 



The equation (fT9~|) is called the gap equation. The solution of the gap equation shows 
the presence of stationary points of the effective action T[(p c }. Thus, the phase structure 
of the theory maybe found by observing the stationary points of r[0 c ] (for a review of 
the quantum field theory in curved spacetime, see Ref. 0). 

3. Scalar Green function at the weak curvature limit 



Following the procedure developed in Refs. [HI El, we solve the Klein-Gordon equation 
and calculate the Green function Gd (|llj) at the weak curvature limit. Note that in 
Green function calculation we actually work in the expansion on the curvatures what is 
known to be appropriate for the study of quantum effects in the inflationary universe 
and also in late-time, dark energy universe. 

The Green function Gd satisfies the modified Klein-Gordon equation, 



-g(x) 



R \ a 



A 



( R 



a-l 



M 2 



iGd(x,x') = 5 (x — x'). 



(21) 



For a practical calculations it is more convenient to introduce G(x, x') and H x which 
are defined by 



G(x, x') 



'(,■) f ^X' 2 iG d {x,x'){-g) l /\x') ( R{X ' ] 



and 



M 2 
R(x^ a/2 



M 2 
( i? V 1 R 



g(.x) 



M 2 



Ar 



a/2 



(22) 



+ a V M 2 



M 2 



X 



R(x) 
~M 2 



-a/2 



In terms of G(x, x') and H x the Klein-Gordon equation ([2"T]) reads 
H x G(x,x') = 5 D {x-x'). 



(23) 



(24) 
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After some calculations the Hamiltonian density H x is rewritten as 

R, 



H x = d^ v d v + a^d» + V x , 
where the potential V x is given by 



(25) 



A 



V x = [(-g)y%(-g)^ dl/ (-g)-^] + £ R-fi+^ 

+ «f [(-.)- 1/4 ^(-.) 1/4 ] - f (f - l) ^ - f ^- 



(26) 



Here we adopt the Riemann normal coordinate expansion ^Hj ■ We set the origin of 
the coordinates at x and write y = x — x' and neglect terms proportional to y a , a > 2. 
Hence, the Eq. ()24)1 is rewritten as 

{V^r + (pd» + + x 2 + PX) G(y) = 5 D (y), (27) 



where 



Cm 

x 2 



a- 



R;fi 

~R 



a: 



R 



R 



R;[iR;v 



2 i ^0 , 2 



1 

" 6 



R 



a /a \ R^R-n a OR 
2\2~) ~ R 2 2~R~ 



= -f4 + u, 



j/=0 



Performing the Fourier transformation 

al D p 



G(y) 



(28) 
(29) 



y=o 
(30) 

(31) 



(32) 



(2tt) d 

one finds that the equation (|27j) becomes 

(V + <^ + + X 2 - 0^) G(p) = 1. (33) 

First we consider the solution for a positive \ 2 . In this case the following proper- 
time form for the solution of Eq. (J3*3*j) is assumed 

( J iSe -X 2 s e P^A^{s)p-+iB^s)p^+C{s) _ 



G(p) 



(34) 



Substituting this into the Eq.(j2HJ) one obtains 

/>oo 

/ als (-p 2 + <X + + x 2 ~ i^X») 

Jo 



X e ~^ 2 s e P^A„ v {s)p v +iB„(s)p^+C{s) 



where 



(35) 
(36) 
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On the other hand, there is an identity 



OO 1 

ds 

ds (Y 2 - rf^p" - z^V - ^ 



ds 



ds 



ds 



x e -x 2s e P fiA ^( s )p" +iB ^ s )p^+ c ( s ) 



= 1, (37) 
Comparing the Eq. (j35|) with the Eq. ()37|) . the simultaneous differential equations are 
obtained 



dA 



ds 
dB„ 



1 — ^p(A p u + A u p ), 
-C-£ w 5 p -/3 P (^ + A/), 



ds 
dC 

~ds~ ~ 

The solutions of these differential equations are found to be 



(38) 
(39) 
(40) 

(41) 
(42) 



B — — T'C + PC 2 + C 2 e^ s - 2PC 1 Cie-*\ 
C = (3r\s - (f3C l ) 2 s + PC^e* - {^fC^ 

+ C 3 , (43) 

where C%, C 2 and C3 are constant matrix parameters which are determined by the 
boundary conditions. In the flat space time, R — > 0, Eq. (j33j) must reproduce a Klein- 
Gordon equation in the Minkowski space. Thus, the functions A, B and C must satisfy 



A -> s, 

B -> 0, 

c ->0, 

at the limit £, £, /3 — > 0. 

These boundary conditions fix the parameters Ci, C 2 and C3 to be 
1 



-A 



-1 



c 2 = -2/?r 2 +r 1 c, 
c , 3 = ^(/9r 1 )(/5r 2 )-/5r 2 c. 

Therefore the functions A, B and C are simplified to 

s= -rta-e-^+^r^i-e^) 2 , 
c= -/?r 2 c(i-e" es -^) 



+ 2(/?r 1 )(/?r 2 ) 



1 



(l- e -^-^)--(l-e- 



2£s 



2^) 



(44) 
(45) 
(46) 



(47) 
(48) 
(49) 

(50) 
(51) 

(52) 
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Since the parameters (, £ and (3 are small enough at the weak curvature limit, we 
neglect the higher order term about these parameters for a small s and get 

A~s-£s\ (53) 
B~-( s + {3s 2 , (54) 

(55) 



C ~ -PCs 2 - -P 2 s 3 ~ 0. 
2 3 



Inserting these solutions into Eq. (j34j) . the Green function G(p) follows 

»A 

G(p) = I dse~ (pli ^" pl '- i t 3 ^) s2 +(p 2 -x 2 -iC^)s ^ 



(56) 



where a scale A ~ \J\R\/M 2 is introduced. It corresponds to the upper limit where we 
can neglect the higher order term about the parameters (, £ and (3. Since the integrand 
becomes smaller and smaller above the scale A, we drop the contribution from the region 
above A. It is straightforward to extend above analysis to a negative \ 2 an d obtain 

-A 

c [ se -(p f *^p v - i i 3 vP f *) s2 +(p 2 -x 2 -iCvP' J ') s ^ (57) 



G(p) = 

Integrating over s in Eq.flHJ) and (JEU), it follows 

(p2 -x 2 - iC^f 



G{p) 



7T 



X 



1 + erf (-sgn( X 2 ) 



cxp 



2 , p 2 -x 2 - 



(58) 



where erf(z) is the error function. In the configuration space the Green function G(y) 
reads 



G{y) 



sgn(x 2 )A 



dse x 



d D p c ipHy U( j>> l A uv p-'+iB u pi' 



At the origin it reduces to 

f>sgn(x 2 )A 

G(0) 



dse x 

sgn(x 2 )A 



(2vr) D 



(59) 



r/.s ( det — 

7T 



-1/2 



-V~W " • (60) 

To perform the momentum integral the Wick rotation maybe applied. Neglecting the 
higher order term about (, £ and (3, we obtain 



G(0) 



(47T 



,D/2 



sgn( X 2 )A rfs 



1 • 



(61) 



Hence, the scalar Green function is found at small curvature. It will be used in the 
calculation of the one-loop effective action. 
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4. The one-loop effective action 

The effective action of the present model is given by (fTTj) . Here we normalize the effective 
action to satisfy 

r[0 = O] = S[<j> = 0]. (62) 
In this normalization the one-loop effective action is given by 

m = S[<f>] - y [In {Det iG^tf)} - In {Det iGf{<? = 0)}] . (63) 
The right hand side of Eq. ()63j) contains a term proportional to 
In {Det iG^i^ 2 )} - In {Det iG d l (4> 2 = 0)} 
= Tr {lniG^ 2 )} - Tr {in zG^(0 2 = 0)} 



Ao 
2 



d v x(-g) 



R 

M2 



dm iGd(x,x;(J) = m). 



(64) 



From the Eq. (j6Ij) the Green function at the weak curvature limit is found to be 



G(x, x; 4> 2 = m 2 ) 



S gn(x 2 )A ds 



*D/2 



1 + ~tr£s ) = m >. (65) 



The relationship between Gd and G is given by 



iGd(x, x; <p 2 = m 2 ) 



1 



R{x) 



The effective action at the weak curvature limit reduces to 

h\ 



G(x, x; 4> 2 = m 2 ). 



Integrating over m 2 one gets 



rf 



4 ( 4vr )D/2 

x / d D x\f^g 



h 



(66) 



(67) 



S[<f>} + 



sgn(x 2 )A ^ s 



dm 



sgn(x 2 )A ds 



,D/2 



1 + -tr£s ) e 



-x 2 (<t> 2 =™ 2 ) 



i + it*.) ( 



s d/2+i\- , r * , • (6* 

Performing the integration about s and neglecting the higher order term about 1/A, we 
get 

r[0] = S[<J>} 

Dtr£ 



+ 



h 



~ ) / d D xV^ 



ix\<P 2 = o)) D / 2 1 



(x 2 (<P 2 )) D/2 i 



£tr£ 

4x 2 (0 2 = 0) 



4x 



The effective Lagrangian density is defined by 



d D x\f^gC 



eff- 



(69) 



(70) 
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From Eqs.((TJ) and (|69p it is found to be 



with 



4p° 



'0 + "y Vfj 



+ 



h 



-r 



1 



2(4tt) d /2 

( Y 2 U 2 = 0) )D/2 ( l ^trf \ 



Dtr£ 



4x 



tr£ 



a 



i2 



i? 2 



10 



aUR 
2~R~ 



(71) 



(72) 



,(73) 



In the limit a — > and i? — * the effective Lagrangian density (J7TJ) reproduces the one 
in the Minkowski space. 



5. The one-loop renormalization 



The effective action ()69j) and the effective Lagrangian density ()71|) are divergent in 
two and four dimensions. To obtain the finite value we must renormalize the theory. 
(Remind that theory is not multiplicatively renormalizable in curved spacetime [5 ). At 
the present order we impose the following renormalization conditions 



#T[0] 



<9 4 r[0] 



4>=Mx 



(j>=M 2 



(- 



(74) 



(75) 



where Mi and Mi are the renormalization scales." 1 " From these conditions one obtains 
the relationship between the renormalized parameters fi r , £ r , A r and the bare ones /xq, 
£o 5 Ao- 



Z r R = f 4-Z R-^M? 



h 



R 



2(4tt) d /2 \M 2 

x {r (i - tf{Ml)) D '^ 



Ar 



l + A M 1 2 (|-2)(x 2 (M 1 2 ))- 1 



trf 



+ For a higher order calculation we must impose the other conditions to renormalize the wave function 
for the field <j> and the cosmological constant. 
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h 



2(4vr) D / 2 \M 2 



(£5) (76) 



A r — Aq 



h 



R 



\ 2 



r ( 2 - 1 ) ( * 2(M2 ))D/2_2 [ 3 + 6A ° M22 (? - 2 ) 



+A 2 M 4 



-3 (X 2 M)- 



+ 2 r(3 



(X 2 (M|)) D / 2 - 3 
D 



3 + 6X MU^-3)(x 2 (M^r 1 



+A 2 M 4 



D 
~2 

4)(X 2 (MI))- 



A + 



2(4vr)^/ 2 VM 2 / "".v— (77) 

By using these renormalized parameters the effective Lagrangian density 1)71)) is 
rewritten as 



eff 



2k 2 



R + 



(- 

1 M 2 



1 



A, 
4! 



+ 



ft 

4(4tt) d / 2 
ft 

2(4tt) d / 2 



Kf{D,M 1 



+ 



ft 



48(4tt) d / 2 



(xW)) D/2 1 



X 2 r g(D,M 2 )(^-6Mt 



2 J,2\ 



4 X 2 (0 2 ) 



(x 2 (o)) D/2 1 



£tr£ 
4X 2 (0)/J 

At the two-dimensional limit Eq.(|78jl becomes 



(7f 



C 



2D 



eff 2k 2 



^0 2 + ^0 2 + ^m 2 2 - 
2 2 4 1 4! 



f 8tt 



2 rV ^ 



x 2 (m 2 ); 



+ X 2 (0)--trC In 



X 2 (0) 



- ^A 2 2 M 2 1 



1 



X 2 (M 2 ) 



~ Aii 



- 6M 2 , 



X 2 (M|) 



3 - 6A r M 



ri "2 2 



X 2 (M|) 



2A 2 M 4 



1 - A r M 2 



1 



X 2 (M 2 ) 



(79) 



r 2 (X 2 (M 2 )) 2 



2x 2 (M|) 



1 - 4A.M 



r 2 2 



X 2 (M|) 

Taking the four-dimensional limit of Eq. (J?8*|) . we obtain 



+ 2A"M; 



2 M 4 - 

r 2 (X 2 (M 2 )) 2 



trf 



^Sf — 7 nR + f 77? 



' e// 2k 2 



\M 2 



'IMS 



r ,4 
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+ 



h 



1287T 2 



A, 



2i2\ 



A r M| 



1 A 2 M 2 4 



X 2 (M|) 6 X 2 (M| 



+ 3A r </>Y(0) + - 2X r 4 2 X 2 {Ml) - 2( X 2 (0)) 2 ln X2(02) 



X 2 (0) 



- 2A r 2 x 2 (O) In 
+ 3A 2 M 2 2 In 
-j|(0 4 -6M 2 2 ) 



X 2 (0 2 ) 

X 2 (M 2 ) 
X 2 (M|) 



2"-- "x 2 (M 2 ) 



2X r M 2 



A 2 M 4 



X 2 (M|) (x 2 (M 2 ))^ 



N 2 A 2 M 2 2 X 2 (0 2 ) 

+ 9/a/r9N - 2% (0) In 0/ns - X r (j) z In 



3(x 2 (M 2 ))3 

xW, 



(80) 



trf 



X 2 (M 2 ) V 7 X 2 (0) r " X 2 (M 2 ) 

Therefore the ultra-violet divergences in the effective Lagrangian density (fTT|) and 
the effective action (J69)) disappear after the one-loop renormalization. Note that to 
make sense to divergent expressions one can apply any other regularization (say, zeta- 
regularization ^1]). Of course, the question of dependence from the regularization in 
such effective models remains. 



6. Phase Structure in FRW spacetime 



It is expected that the non-linear curvature coupling in our model leads to non-trivial 
consequences for the phase structure in a non-static spacetime. (This may have the 
interesting applications for early-time, accelerating and late-time, accelerating universe 
when quantum gravity effects are negligible). Here we consider the model in the D- 
dimensional Fridmann-Robertson- Walker spacetime. By spirit it is similar to the study 
of phase structure for scalar self-interacting theory, scalar electrodynamics and GUTs in 
curved spacetime (see [121 EH EH EEl EG3 El EE] and for complete list of refs. and review 
j^]) where curvature-induced phase transitions were discovered. The FRW spacetime is 
defined by the metric 

dr 2 



ds 2 = dt 2 - a 2 (t) 



+ r 2 dn D _ 2 



1-Kr 2 

where K corresponds to the curvature of the D 
R in this spacetime is given by 



(81) 



1-dimensional space. The curvature 



R = 2(D-1) 



D 9 D 
H + — H 2 + — 
2 



2K 



\2) 



where the Hubble rate H is defined by H 



R ]0 = 2{D - 1) 



H + DHH — (D — 2)H 



a/ a. The first derivative of the curvature is 
K 



(83) 



R. t = 0. 



(84) 
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Thus, it follows 



R^R ;fl = A(D - l) 2 



H + DHH — (D — 2)H— 

a 2 



DR =2(D-1) H+(2D-l)HH + D(D-l)HH 2 + DH 2 



(D — 2) (dH 2 + H-3^\ 



H\ K 

2 



55) 



56) 



Here the spacetime with if = and a(t) = a t ho is considered. In this spacetime 
the curvature invariants are 

R =2(D-l)h ^jh -lJ^, (87) 

= 1 

R 2 t 2 ' 

^ = -2[(D-l)/io-3]i. (89) 

The ground state of the theory depends on the spacetime structure and the 
parameters of the theory under consideration. To study the phase structure of the 
<f) 4 theory in FRW spacetime we evaluate the effective action, C e ff. The expectation 
value of the field is obtained by searching of the stationary points of the effective 
action. The effective action develops a small imaginary part for a negative x(0)- We 
evaluate a real part of the effective action to find the ground state. 

In an expanding universe the scale factor a(t) increases as the time runs. Here 
we suppose that the scale factor a(t) is proportional to the time t, i.e. a(t) = at 
and numerically calculate the renormalized effective action (|78|). All mass scales are 
normalized by an arbitrary mass scale M and h — 1. The renormalization scale is 
chosen to be M 1 = and M 2 = O.IM^" 2 '/ 2 . 

First we consider the stationary and spatially homogeneous <fi. In this case the 
kinetic term of </> disappears. One can define the ordinary effective potential V(4>) by 

V{<j>) = -C eff . (90) 

The expectation value of the field <fi is defined by the minimum of the effective potential. 

We illustrate the behavior of the effective potential for a minimal gravitational 

coupling £ = at /i r = 0.1M in Figs.l, 2 and 3. Here we set D = 3.8 to regularize 

the theory with keeping the general covariance. It corresponds to the cut-off scale, 

A ~ 500M 2 , where M 2 is a renormalization scale. Since the spacetime curvature behaves 

as t~ 2 , it raises up for a small t. It is not valid to apply the weak curvature expansion 

at the beginning of the universe. In Figs. 1, 2 and 3 we solve the gap equation only 

for 2 < tM. As is shown in the figures the effective potential seems to be almost 

proportional to R a oc t~ 2a . 

In Figs. 1 and 2 the expectation value (<j)) is located near the stationary point of 

6 

the classical potential, <p 2 = —fi 2 . For the minimal gravitational coupling the classical 

A 

part of the effective action has no curvature dependence except for the overall factor 



The one-loop effective action in 4 theory coupled non-linearly with curvature 



14 



0.0002 


-0.0002 
^ -0.0004 
^-0.0006 
of -0.0008 
°- -0.001 
-0.0012 
-0.0014 
-0.0016. 



t=20/M 
t=30/M 
t=40/M 
t=50/M 



0.2 0.4 0.6 0.8 

cp/M( D - 2 ) /2 




40 60 
tM 



100 



Figure 1. Behaviour of the effective potential and the mass gap for a = 0, ho = 1, 
£ = 0, fi r = 0.1M, A = 0.1M 4 -° and D = 3.8. 
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Figure 2. Behaviour of the effective potential and the mass gap for a 
£ = 0, fi r = 0.1M, A = 0.1M 4 - D and L> = 3.8. 
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Figure 3. Behaviour of the effective potential and the mass gap for a — 1, ho = 1, 
£ = 0, = 0.1M, A = 0.1M 4 - D and D = 3.8. 



The one-loop effective action in c/> 4 theory coupled non-linearly with curvature 



15 



3x10 39 
^ 2x10 39 
> 1x10 39 

CD 

= 
-1x1 39 
-2x10 39 




t=4.0x109/M 
t=4.5x109/M 



1x10 

cp/M< D -2) /2 



3.5x1 9 



5.5x1 9 
tM 



7.5x1 9 



Figure 4. Behaviour of the effective potential and the mass gap for a = 0, ho = 1, 
£ = 0, fi r = 0, A = 0.1M 4 " D and D = 3.8. 
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Figure 5. Behaviour of the effective potential and the mass gap for a = 0, ho = 1, 
£ = Cconformal, Mr = 0, A = 0.1A/ 4 - D and D = 3.8. 



and the Einstein-Hilbert term. On the other side, the radiative correction depends on 
the curvature. At a = our model reduces to the ordinary c/> 4 theory. As is known, the 
radiative correction has only a small contribution to spontaneous symmetry breaking 
for a large fi r in the ordinary theory. In the case of a negative a the classical part of 
the effective action (fT%|) is enhanced as the curvature decreases. Contributions from 
the radiative corrections are lost at a large t. For a positive a the classical part of 
the effective action is suppressed except for the Einstein-Hilbert term as the curvature 
decreases. Hence we observe the curvature dependence of the expectation value ((f)) only 
if a > 0, see Fig. 3. 

Next we plot the behaviour of the effective potential for /i r = in Figs. 4, 5, 
6 and 7. Here we consider the minimal and conformal gravitational coupling, £ = 
and £ = Conformal = (D — 2)/[4(_D — 1)] respectively. (Of course, the scalar sector is 
conformally invariant for correspondent choice of £ only when a = 0.) For \i r = 
spontaneous symmetry breaking can not take place on the classical level of the ordinary 
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Figure 6. Behaviour of the effective potential and the mass gap for a = 1, ho = 1, 
£ = 0, Mr = 0, A = 0.1M 4 ~ D and D = 3.8. 




(p/M(D-2)/2 tM 
Figure 7. Behaviour of the effective potential and the mass gap for a = 1, ho = 1, 

£ - Cconformal, Mr = 0, A = 0.1A/ 4 - D and D = 3.8. 

(f) A theory. The radiative correction has an essential role to break the symmetry in the 
ordinary theory, i.e. a = 0. It is known as a radiative symmetry breaking. First order 
phase transition is observed in Figs. 4 and 5. However, the radiative correction has only 
a small effect. The expectation value generated by the radiative symmetry breaking is 
extremely small, as is shown in Figs. 4 and 5 at a = 0. A long time is necessary until 
phase transition takes place. 

In the case of a positive a the ratio of the radiative correction to the classical part 
of the c/> 4 Lagrangian is enhanced as the curvature decreases. In Figs. 6 and 7 we 
also observe the first order phase transition and find that the expectation value (0) is 
extremely enhanced in comparison with the case a = 0. For a = — 1 only a symmetric 
phase is realized at /i r = and then radiative symmetry breaking does not occur. 

As is shown in Figs. 6 and 7, the expectation value of <j) has non-negligible time 
dependence for the case \i = 0. Hence we next consider the spacialy homogeneous but 
time dependent at \i = and a = 1. Then we should solve the gap equation (JTHJ). 
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From Eq. (|19|) the equation of motion for is obtained 

1 8 f lD , „ dVU) 1 n . (R\ a „u, 



dV{4>) _( R 
~dd) VM 2 



0. (91) 



To find an exact solution one needs to solve this equation for a general time-dependent 
form of 0. However, it is instructive to consider the solution of Eq. (j91j) for a special 
form of 0. In the present paper it is assumed that 

2 

0(t) = (0(t)> = vt\ = |r<0W> 2 , (92) 

where v is a constant parameter. Here we fix the parameter x at 1/2 or 1 and numerically 
solve the equation of motion for 0. 

Under the assumption (|9*2"|) Eq.(|9*T|) reads 

[(x - 1) + (D - l)h - 2a] + (J^j = o. (93) 

Thus the expectation value (0) is found by observing the stationary point of the following 
function 

F[0] = [(x -1) + (D- l)h - 2a] ^(tf + (JL^j " V {<j>). (94) 

The solution is shown in Figs. 8 and 9. As is seen in Figs. 8 and 9, we observe the 
first order phase transition again. In the case of the minimal gravitational coupling 
we find the two steps of the first order phase transition for x = 1. In Fig. 8 (a) the 
solution becomes unstable for tM < 15.5. The weak curvature approximation and/or 
the assumption (4>{t)) = vt can not be applicable for such a small t. In all the cases the 
mass scale v depends on the time t. However, we observe that the mass scale v is almost 
static between tM = 80 and tM = 90 for x = 1 and between tM = 90 and tM = 100 
for x = 1/2. It is expected that there is a solution with gradually decreasing x after 
the first order phase transition near the time tM = 60. On the other hand, we observe 
in Figs. 6, 7, 8 and 9 that the time of the first order transition near tM = 60 has no 
large dependence on x in the present analysis at a = 1 and h = 1. It seems to be a 
characteristic feature of the theory with a = 1 in the FRW spacetime. 

Hence, in the present section we numerically estimated phase structure for scalar 
theory in curved spacetime under consideration. The phase structure strongly depends 
on the sign of a. Compared with the usual 4 theory in weakly curved space, i.e. 
a = 0, a positive a significantly enhances the symmetry breaking. For a negative a the 
radiative symmetry breaking does not take place. As it was shown, mainly first-order 
curvature induced phase transitions may occur here. 

In all figures we put D = 3.8 and ho = 1. The spacetime curvature is given 
by Eq. ()S7jl and found to be R ~ 5.0/t 2 . The typical scale found to be tM ~ 50, 
(0) ~ 0.5M( D - 2 )/ 2 in Figs. 1-3, 6-9 and tM ~ 10 10 , (0) ~ 10- 9 M( D " 2 )/ 2 in Figs. 4, 
5. The curvature scales for such time scales are given by R ~ 0.002M 2 and 10 _20 M 2 
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Figure 8. Behaviour of the mass scale v for a = 1, ho = 1, £ = 0, fj, r = 0, 
A = 0.1M 4 ~ D and D = 3.8. 
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Figure 9. Behaviour of the mass scale v for a = 1, ho = 1, £ = Cconformab Mr = 0; 
A = 0.1M 4 - D and D = 3.8. 



respectively. If we consider the electroweak scale Higgs boson, ((f)) ~ (lOOGeV)^^" 2 ^ 2 , 
we obtain the typical mass scale M ~ 200GeV in the former case and M ~ 10 11 GeV in 
the latter case. The cut-off scales which corresponds to D = 3.8 are A ~ ITeV, near the 
scale of a large extra dimension, and 10 12 GeV, respectively. For both cases the typical 
scale of time and curvature are given by t ~ 0.1(GeV _1 ) and R ~ 100(GeV 2 ). Thus the 
quantum correction induced by the curvature effect has an interesting contribution to 
the effective action at the very early stage of the universe. The time scale of the phase 
transition depends on the coupling constant A. For a stronger coupling the radiative 
correction is enhanced and a smaller time scale is obtained. 

Thus we consider the regularized theory and calculate the effective potential in a 
D-dimensional spacetime with D = 3.8. As is shown in the previous section, we obtain 
the finite effective potential (JHUj) in four spacetime dimensions. Here we show some 
results in four dimensions. Evaluating the effective potential (|8T)|) we draw the solution 
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Figure 10. Behaviour of the effective potential and the mass gap for a = 1, ho = 1, 
fir = 0, A = 0.1 in D = 4. 



of the gap equation. The first order phase transitions is observed as is shown in Fig. 
10. Compared with the case D = 3.8 (Figs. 6 and 7) the critical time becomes late and 
the the two steps transition does not realised in four dimensions. 

The main lesson of this consideration is that it is unlikely that radiative corrections 
may strongly influence the dark energy universe in such theory. However, quantum 
effects may become quite important for phantom universe where Big Rip may occur 
[T§1 I2U] . It was already demonstrated that in such circumstances already quantum 
effects of free fields may prevent Big Rip occurrence [2U 1221 1231 121] • However, in order 
to analyse this problem in our framework one needs the calculation of the one-loop 
effective action at strong curvature. 



7. Dynamical mechanism to solve the cosmological constant problem 



7.1. Dynamical cosmological constant theory: exact example 

There are some proposals to solve the cosmological constant problem dynamically [2]. 
One of the possible solutions of the cosmological constant problem is pointed out by 
Mukohyama and Randall In [3], the following type of the action similar to the one 
under consideration has been proposed: 



d^xsj—g 



R_ 2 (K%<p&><p) q 

2k 2 i- i- 2qK A f(R)^- 1 



V{<p) 



(95) 



where f(R) is a proper function of the scalar curvature R. When the curvature is small, 
it is assumed f(R) behaves as 

f(R) ~ {Jtf) m . (96) 

Here m is positive. When the curvature is small, the vacuum energy, and therefore the 
value of the potential becomes small. Then one may assume, for the small curvature, 
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V(<p) behaves as 

V(<p) ~V ((p- cp c ) . (97) 

Here Vq and ip c are constants. If q > 1/2, the factor in front of the kinetic term 
of (p in (f9"K|) becomes large. This makes the time development of the scalar field cp 
very slow and it is expected that (p does not reach (p c . This model may explain the 
acceleration of the present universe. The model ()95)1 is also expected to be stable under 
the radiative corrections. In fact, when the curvature is small and the time development 
of the curvature can be neglected, if we rescale the scalar field ip as ip -> RTto-Qhip, 
the curvature in the kinetic term can be absorbed into the redefinition of ip and there 
appear factors including R m ^ 2 i~ l ) i 1 . Therefore if m{2q — l)/q > 0, the interactions 
could be suppressed when the curvature is small and there will not appear the radiative 
correction to the vacuum energy except the one loop corrections. 

We now show that there is an exactly solvable model, which realizes the scenario 
in [Sj. Assume m = q = 1 and (|9T)j) and (J9*7)) are exact: 

f(R) = (3R 2 , V(<p) =V ((p- <p c ) . (98) 

Here /3 is a constant. R 2 term is neglected by putting a = in (|93j) since the curvature 
is small. Then by the variation over the metric, one obtains 



0= - 



+ 



1 f R d p <pd p (p 



-V ((p- (p c ) > Qyn, + 



Rav d p <pd p ipR 



2 I 2k 2 2(3R 2 



2k 2 



f3R 3 



2(3R 2 



On the other hand, the variation over ip gives 







-d„ 



-9 



(3R< 



We assume the four dimensional FRW metric with flat spatial part: 
ols 2 = alt 2 - a(t) 2 ( dxi T ■ 



(99) 



(100) 



(101) 



=1,2,3 



It is chosen that (p only depends on the time coordinate t. Then the (t, /^-component of 
(HE! gives 



3H 2 



144/3 [H + H 2 



+ V (<p- <p c ) 



+ H 



dt 



12/3 [H + H 2 



H + H 2 



J 



72(3 [H + H 2 



3 ' 



(102) 



where H is the Hubble rate. On the other hand, (jlOOj) gives 



/ 



= - 



dt 



a 3 (p 



\ 



V 



36/3 [H + H 2 



-V a 3 



(103) 
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A solution of Eqs. (|1U2|) and (|103|) is given by 

a = a t h « = y) , P = Vc + |r > ( 104 ) 

when ho > 0. In case h < 0, this scale factor a = aot ho f|104|) does not express 
expanding universe but shrinking one. If the direction of time is changed as t — > —t, the 
expanding universe emerges with scale factor a = ao (-t) ho . In the expression, however, 
since h is not always an integer, t should be negative so that the scale factor should 
be real. To avoid this apparent difficulty, we may further shift the origin of the time 
as t — > — t — > t s — t. Then the time t can be positive as long as t < t s . When ho < 0, 
instead of (|1U4|) . one may propose 

a = a (t s -t) h0 (h = -^-) , v = Vc+—^-i, (105) 

The assumption (jlU4j) or (jlU5j) reduces Eqs. (jl02|) and (jl(J3j) to the algebraic equations: 

3/i| <p 2 h 
2k 2 36/3 {-h + 2h 2 f 
cpo (3hp + 1) 
18(3 (-ho + 2h 2 y 



= - ^ + ^ 5 + V <p , (106) 

°-^ 2 36/3(-/i + 2/i 2 ) 3 1 J 

0= ^:°} 3h ° + !l, 2 -Vo, (107) 



which gives 



3 A 4 



54/3 (-1 + 2fe ) fei 
(12h 2 - 2/i - 1) 



3h + 1 , , 

Vo = ± ° . 108 

v/6k 2 (12h 2 - 2h - 1) (-1 + 2ho) 

Since <£>q should be positive, one finds 



u a n I" V13 , 1 + V13 1 
when /? > , — — — < h < — — — or ho > - 



1 — v 13 1 + v 13 I 

when < , h < — or — < h < - . (109) 

12 12 2 

By properly choosing the parameter /3 and Vo, we can always obtain a solution as in 
(gDD or ifTHSjl . For example, if 

». = -i, (HO) 

which gives the following value of the parameter w of the equation of state 

w = -l + ^- = -1.025 , (111) 

we find 

19 /l~5" 

«V = ±— \ — = ±0.388722... . (112) 
34 V 31 v ' 

It is interesting that the value of w (jlll|) is consistent with the observed one. 
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For ho > case, corresponding to (|1U4|) . since R = 6H + 12H 2 , the curvature R 
decreases as t~ 2 with time t. Eq. (jl(J4|) tells that <p approaches to ip c but does not arrive 
at ip c in a finite time, as expected in [Hj. 

As seen from the expression of H in (I1U4|) or (|105|) . if we substitute the value of the 
age of the present universe 10 10 years~ (10 _33 eV) _1 into t or t s — t, the observed value of 
H could be reproduced, which could explain the smallness of the effective cosmological 
constant A ~ H 2 . Note that even if there is no potential term, that is, Vq = 0, when 
(3 < 0, there is a solution 

1 1 - Vl3 . , 

h = -3 < J = -0.2171... , (113) 

which gives the parameter w of the equation of state: w = —3, although w is not 
realistic. 



7.2. Cosmological constant problem in scalar theory non-linearly coupled with curvature 



Here we consider the solution of the cosmological constant problem in the theory 
As is shown in the previous section, one could expect that the radiative correction is 
suppressed when the curvature is small if a < 0. As in the case of if we rescale by 
— > (R/M 2 )~ a / 2 (j), the kinetic term becomes standardly normalized one and does not 
include the curvature. The 4 interaction terms, however, include the factor (R/M 2 )~ a 
and may be suppressed for the small curvature. Hence, there would not appear the 
radiative correction to the vacuum energy except the one-loop corrections. 

As an example of the solvable case, we put a = — 1 and /ijj = 0. The variation over 
the metric gives 



0= - 



1 

2 

M 2 
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2^2 



~R~ 
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A 



9iw + 
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[IV 



2k 2 



+ 2R 9 ^ 



M 2 R 
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Scalar field equation is 







1 



--d„ 



-9 



-9- 



~R 



A 
3LR ( 



(114) 



(115) 



By assuming the FRW metric (jlOljl and 
forms: 



4>(t), Eqs. ()114j) and ()115jl have the following 







3H 2 M " (l^ 2 - ¥') _M 2 (H + H 2 ) ( 



+ b 
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A <4 



(h + 2H 2 J 



(116) 
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a 3 4> 
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dt\H + 2H 2 
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3! [H + 2H 2 



+ £0a 3 . 



With the Anzats 



H 



when hn > or 



t 



H 



h 



U-t 



U-t 



when ho < 0. Eqs. ()116|) and ()117|) reduce to the algebraic equations: 

M 2 fl .„ A 



0= - 



3/ig 
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2k 2 ' 12 (-/i + 2/ig) 



4!^ 



M 2 
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12 (-/i + 2h 2 Y 
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which give 



A 



2 (-/i + 2h 2 ) 36 (-/i + h 2 ) 
3 

= /:2 / 8-9/ip (4-7feo)C 

K {2A(-h + h 2 ) 2 (-h + 2h 2 )h 
-- -6K 2 h {l-2(l-2h )£} 

8-9h (4-7fc K 



(117) 



:il8) 



(119) 



(120) 
(121) 

(122) 



(123) 

2A(-h + h 2 ) 2 {-h + 2h 2 )h j 

Therefore, with the proper choice of parameters, the solution in the form of ()118|) or 
(I119J) follows. We should note the behavior of H f!118|) or ()119j) is almost identical with 
that in (ITU41 or (fTUoT) . 

One may consider the model including two scalar fields ip and <fi as 



d A x\/—g 



2k 2 2(3R 2 



M 2 



2 M 



-V ((p- <p c 
A 



(124) 



Even in this case, for FRW metric ()101|) and assuming ()104J) or (jl(J5j) and/or (|118|) or 
(|119|) . we obtain ()117|) and (JUT)) . The equation corresponding to ()116j) and/or ()120|) 
has the following form: 
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M A 2 A 4 



+ U(-h + 2h 2 ) \2 yo W° 

M 2 (-h + h 2 ) (1, 2 A A M 2 £ , 



12(-/ l0 + 2/ iS rV2"° + 5^ + 4 «■ (125 » 
Then by properly choosing the parameters, we may obtain an exact solution for 
cosmological constant again. 

Hence, we demonstrated that the mechanism proposed in ref.[3] to dynamically 
solve the cosmological constant problem is naturally realized also for the class of models 
investigated in this work. 

The final remark is in order. In [2H] (see also |2fij). the cosmological constant 
damping in the model contained <p 2 R term in the Lagrangian has been investigated. It 
has been shown the possibility to solve the cosmological constant problem. Although 
the problem could be solved in the model, the effective Newton constant G e fr suffers 
the correction like G e s ~ Go/ (1 + a 4> 2 ) with a constant a. Here Go is the bare Newton 
constant. In the solution of model [22], 4> behaves as a linear function of the time t. 
Naturally t could be the age of the present universe, t ~ (10~ 33 eV) , the correction 
becomes very large and the effective Newton constant becomes unnaturally small. In 
the model under consideration (JIJ) with (J2J), similar term is included but due to the 
factor R a , if a > 0, the correction is suppressed since R ~ (10- 33 eV) and R is very 
small. In general a can be also negative. In the solvable case in ()114|) . the correction 
could be R- 2 {\d p <pd p <i) + ±<p A }. Since the R ~ t~ 2 or (t„ - t)~ 2 , by combining (ITTTfl) 
or (|119j) . the correction takes the order of unity. Then the effective Newton constant 
could receive the correction but it is finite, which does not lead to the instabilities of 
the above sort (smallness of gravitational constant). 



8. Conclusion 



In summary, we studied the scalar self-interacting theory whose Lagrangian non-linearly 
couples with some power of the curvature as the effective field theory. The one-loop 
effective action at weak curvature is found here using specific regularization scheme. 
The phase structure of the theory is carefully investigated, the examples of curvature- 
induced phase transitions are presented numerically (for the case of positive power of 
curvature). The comparison with the ordinary scalar self- interacting theory in curved 
space is done. For negative value of a the radiative symmetry breaking does not take 
place. It remains to develop the same one-loop formulation at strong curvature as in 
such a case the phase transitions maybe considered in the vicinity of Big Rip. 

The dynamical mechanism to explain the origin of the extremely small cosmological 
constant in the spirit of papers [21 E] is presented. It is shown that in the class of 
the models under consideration (including the case with two scalars) such a proposal 
maybe realized quite successfully. Unfortunately, despite the preliminary expectations 
the radiative corrections for the theory under consideration are not relevant in the study 
of the dynamical cosmological constant. 
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There are many directions where our approach maybe generalised. In particularly, 
it is known that phase structure of NJL model in curved spacetime is quite rich (for 
a review, see |2Z|)- Hence, it would be extremely interesting to consider the NJL-like 
models non-linearly coupled with power of the curvature in the same spirit. Such models 
maybe even more attractive due to fact that 1/N expansion is more reliable than one- 
loop approximation in non-renormalizable theories. One can also consider Yang-Mills 
theory (or Born-Infeld theory) coupled to curvature powers in the same way. Definitely, 
this may bring the number of the interesting applications in such issues as dark energy 
and cosmo logical constant. 

From another point, the new matter-gravity coupling pQ maybe also considered as 
kind of modification of gravitation itself. In such a situation, when theory somehow 
deviates from General Relativity, it is not clear apriori if the metric formalism and 
first-order formulation should lead to the same physical picture |2H]- (Indeed, unlike 
to General Relativity where metric and Palatini formulations basically coincide, it is 
not the case for our theory (1)PQ| already on classical level [2H])- It would be extremely 
interesting to develop the one-loop effective action formulation in the same way as the 
one developed in this paper when an external gravity is formulated in Palatini form. 
This will be done elsewhere. 
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